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ABSTRACT 

Effects of magnetic field on stellar differential rotation are studied by comparing magnetohydrodynamic 
(MHD) models and their hydrodynamic (HD) counterparts in the broad range of rotation rate and in varying 
initial rotation profile. Fully-compressible MHD simulations of rotating penetrative convection are performed 
in a full-spherical shell geometry. Critical conditions for the transition of the differential rotation between faster 
equator (solar-type) and slower equator (anti-solar type) are explored with focusing on the “Rossby number 
(Ro)" and the “convective Rossby number (Ro con v)"- It is confirmed that the transition is more gradual and 
the critical value for it is higher in the MHD model than the HD model in the view of the Ro con v-dependence. 
The rotation profile shows, as observed in earlier studies, the bistability near the transition in the HD model, 
while it disappears when allowing the growth of magnetic fields except for the model with taking anti-solar 
type solution as the initial condition. We find that the transition occurs at Ro ~ 1 both in the MHD and HD 
models independently of the hysteresis. Not only the critical value, the sharpness of the transition is also 
similar between the two models in the view of the Ro-dependence. The influences of the dynamo-generated 
magnetic field and/or the hysteresis on convective motion are reflected in the Ro. This would be the reason why 
the transition is unified in the view of the Ro-dependence. We finally discuss the Ro-dependence of magnetic 
dynamo activities with emphasis on its possible relation to the kinetic helicity profile. 

Subject headings: convection-magnetohydrodynamics (MHD) - Sun: interior - Stars: rotation 


1. INTRODUCTION 

The differential rotation (DR) is believed to be a key in¬ 
gredient in organizing large-scale magnetic fields in the solar 
interior. Most of standard solar dynamo scenarios rely heav¬ 
ily on, so-called fTeffect as the amplification process of mag¬ 
netic fields to reproduce observed solar magnetic activity with 
cyclic polarity reversal s and butterfly-shaped s patiotemporal 
migrations (see, e.g., Charbonneau 2005, 2010 for reviews). 
However, we have not yet arrived at a full understanding of 
the physical mechanism for maintaining the solar differential 
rotation. See, e.g., |Miesch| ( [2005 [ ) for a review. 

The internal rotation profile of the Sun has been uncovered 
by the global helioseismology (e.g ., [Kosovic hev et al.( l997 ; 
[Thompson et al.|200^|Howe|2009| ). A striking feature of the 
solar rotation profile is its equatorial acceleration: The plasma 
near the equator rotates faster than other parts in the conven¬ 
tion zone. While the bulk of the convection zone is character¬ 
ized by conical iso-rotation contours, the strong radial shear 
is concentrated in the near-surface layer and the tachocline 
underlying the convection zone (e.g.,|Schou et al.|19 98 ; Char- 


bonneau et al.||1999[|Basu & An tia 2001). The equatorward 

angular momentum transport due to the turbulent Reynolds 
stress induced by the rotating stratified convection is mainly 
respo nsible for the solar-type DR with equatorial acceleration 
(e.g.JKrause & Rudigerll9741 |Ruediger|l 9891 |Kitchatinov &| 

|Rudiger|l995 ; Miesch 200$f 

possibility of the opp< 


The possibility of the opposite profile, or anti-solar type 
DR, in which the equatorial region is rotating slower, has 
been suggested numerically since the early 3D hydrody namic 
simulation of rotating spherical-shell convection by [Gilman 


( |1977| ) over a b road range of parameters and various Sim ula 
tion setups (see Gla tzmaier & Gilman 1 982|jBrun & Toomre 

& Frey 


2002[ Aurnou et al.| 2UT)7[ |Steffen & Freytag||2007[ |Chan 
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2010; |Kapyla et al. 2011| Gastine et al. 2013). However, 
although the solar-type DR has been frequently observed in 
various types and ages of the stars, the anti-solar type DR has 
only been reported for a few of K giant stars observed with the 
Doppler imaging technique (Strassme ier et al.|2003 |Weber et 


|al.|2005||Kovari et al.|20 07 1 2014]). It is still unclear how com 

mon the anti-solar type DR is in the solar-like main-sequence 
dwarfs and what really separates the two rotation regimes. 

Bearing in mind these situations, the transition between the 
solar and anti-solar type rotation profiles has been extensively 
studied by the simulation of spherical shell convection in re¬ 
cent years. [Kapyla et al.[ ( [201l| ) found in fully-compressible 
convection simulation that the transition is characterized by 
the Coriolis number (inverse of the Rossby number which is 
defined in §2 ) and is only weakly d ependent on the density 
stratification. Guerrero et al. ( |2013| ) explored the physics of 
large-scale flows in solar-like starsTiy using an anelastic large- 
eddy simulation with stratification resembling the solar inte¬ 
rior. They confirmed that the two regimes of the DR exist even 
in the solar-like strongly stratified internal structure. The crit¬ 
ical Rossby number they obtained was consistent with that in 
Kapyla et al.]( |2011| ) 


The systematic parameter study of the transition between 
the sol ar an d anti-solar type rotation profiles was conducted 
by [Gastine et al.| ( j2Q14 ) with rotating spherical shell Simula 
tions of anelastic and Boussinesq convections (see also |Gas 
tine et al. 2013| ). They found, by combining their massive 
simulation results with the models of the different research 
groups, that the transition is controlled by the “convective 
Rossby number", almost independently of the model setup, 
such as thickness of the convective envelope and density strat¬ 
ification. In addition, they found that two kinds of DR profiles 
are two possible bistable states around the transition, suggest¬ 
ing the hysteresis of the stellar rotati on pr ofile. The bistability 
of the rotation profile, disco vered byfG astine et ak] ( [2Q14[ ), was 
confirmed independently by [Kapyla et al.| ( [2014[ ) in hydrody- 
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namic simulations of the compressible convection. However, 
the influence of magnetic field on the rotation profile and its 
transition is still controversial despite the magnetic field is an 
inevit able outcome of electrically-conduc ting fluid motions 
(e.g. ,|Moffatt| 1978]JKrause & Raedler| 1980| ). 

In|Masada et al.|(J2013J) (hereafter MYK13), we suggested, 
for the first timeTthat the formation of the solar-type DR is 
associated with the development of the magnetic field : It 
was found that the anti-solar type DR profile at the early dy¬ 
namo kinematic stage transits to the solar-type profile after 
the dynamo-generated magnetic field beginning to affect the 
convective motion (see §4.1 in MYK13). Such a magnetic in¬ 
fluence on the rotation profile was confirmed by independent 
groups recently. |Fan & Fang|(j2014| ) performed anelastic con¬ 
vective dynamo simulation of the model with the realistic so¬ 
lar internal structure, solar rotation rate and solar luminosity. 
They showed that the dynamo-generated magnetic field plays 
a crucial role in attaining the solar rotation profile in the ac¬ 
tual solar parameter regime, i.e., the rotation profile becomes 
the anti-solar type without the magnetic field. They also re¬ 
ported that the bistability of the rotation profile, discovered 
in earlier studies, disappears when allowing the evolution of 
the magnetic field. The facilitation of the solar-type DR and 
the disappearance of the bistability of the rotatio n pro file due 
to the magnetic field were confirmed by [Karak et al.| ( |20l4| ), 
in which they studied the effect of the magnetic field oiTthe 
rotation profiles around the transition by fully-compressible 
convective dynamo simulations. 

In this paper, we systematically study the effect of the mag¬ 
netic field on rotating spherical shell convection and resultant 
mean DR by comparing the MHD models with their hydro¬ 
dynamics counterparts (HD models) in the broad range of the 
initial rotation rate. The bistable nature of the rotation profile 
is revisited under the influence of the magnetic field. A self- 
consistent, fully-compressible Yin-Yang MHD code which 
was developed in MYK13 and a stellar model consisting of 
the convection zone and the radiative zone are used for the 
simulation. The primary objective of this paper is to explore a 
key parameter which controls the transition between the anti¬ 
solar type and solar-type DR profiles and its dependence on 
the magnetic field and the initial rotation profile (§3.1-3.3). 
We cast a spotlight on two diagnostic parameters, i.e., the 
“Rossby number" and the “convective Rossby number" (see 
§2 for the definitions). The magnetic dynamo activity in the 
MHD models with the different rotation rates is also studied 
in §3.4. After discussing the relation between the magnetic 
dynamo activity and the convection properties in §4, we sum¬ 
marize our findings in §5. 

2. SIMULATION SETUP 

The simulation setup is almost the same as the model A 
of MYK13. We numerically solve a MHD convection sys¬ 
tem in a spherical shell domain defined by (0.6R < r < R), 
(0 < 6 < 7 r), and (— 7 r < (j) < i r), where r, and are ra¬ 
dius, colatitude, and longitude, respectively. Our model con¬ 
sists of two-layers, which qualitatively resemble the solar in¬ 
terior: stably stratified layer of thickness 0.1 R in the range 
(0.6 R < r < 0.1 R) and surrounding convective envelope of 
thickness 03R in (0.1R <r<R). 

The basic equations are the fully-compressible MHD equa¬ 
tions in the rotating frame with a constant angular velocity 
fl = Qoe z which is parallel to the coordinate axis (0 = 0): 


— =-V-/ 
dt J ’ 


( 1 ) 


£ £=-V-(vf)-Vp+jxB 


+pg + 2pv x f l + v 
dp 

— = -v -Vp-^pS/ •v 


1 


V' v+ -V(V -v) 


dA 

dt 


+( 7 - 1 ) [v-(Kvr)+ 7 ?j 2 +$] , 
= v x B — gj , 


( 2 ) 

( 3 ) 

( 4 ) 


with 


<f> = 2 z/ 


e ij e ij ^ 


_ 1 / dvi dvj 
’ €lj 2 \dxj +dxi 


9 = ~go/r 2 e r , B = V x A , j = V x B . 

Here the mass density p , pressure p , mass flux / = pv , mag¬ 
netic field’s vector potential A are the basic variables, go is 
the gravitational acceleration, assumed to be constant. We as¬ 
sume an ideal gas law p = ( 7 - 1 )pe t with 7 = 5/3, where e* 
is the internal energy. The viscosity, electrical resistivity, and 
thermal conductivity are represented by v, 7 , and k respec¬ 
tively. 

The initial condition is a hydrostatic equilibrium which is 
described by a polytropic temperature distribution with the 
polytropic index m 


dT _ -gp/r 2 
dr c v ( 7 ~l)(m+ 1 ) 


( 5 ) 


where c v is the specific heat at constant volume. We choose 
m = 1 and 3 for the upper convective envelope and the lower 
stable layer, respectively. The thermal conductivity k is de¬ 
termined by requiring a constant luminosity, L, defined by 
L = -47wr 2 dT / dr , throughout the domain. 

Non-dimensional quantities are defined by setting R = go = 
po = Po = 1 where po is the initial density at r = 0.6R. The 
units of length, time, velocity, density, and magnetic field are 

R, y/R/go, yfgoR, Po and yfgoRpoPo, respectively. The defini¬ 
tions of the Prandtl, magnetic Prandtl, and Rayleigh numbers 


are 


v v god 4 

Pr = Pm = Ra = 


X 


T] 




1 ds 
c p dr 


( 6 ) 


where x = K m/( c pp m ) is the thermal diffusivity measured at 
the middle of the convection zone (r = r m ), and d = OAR is 
the thickness of the convective envelope. The stratification is 
controlled by a normalized pressure scale height at the surface 
£0 = < 7 ( 7 - l)T s /(go/R), where T s is the temperature at r = 
R. Here we use in all the models £0 = 0-3, yielding a small 
density contrast between top and bottom boundaries about 3 
(see Figure 1 in MYK13). This is a major difference between 
our model and actual Sun. 

In this paper, we focus on the two physical non-dimensional 
quantities, one is the “Rossby number" (ratio of inertia to 
Coriolis forces) and the other is the “convective Ro ssby n um 


ber" (ratio of buoya ncy to Coriolis forces) (e.g., |Gilman|1977 
Gas tine et al.|2014 ) which are given by 


R 0 = Co -1 = 


'Yms^-7 

2 fi 0 


ROc, 


-(&) 


1/2 


( 7 ) 


respectively, where ki = 2 k/ d is the wavenumber of the largest 
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TABLE 1 

Summary of the simulation runs. The labels of MHD models contain the letter M. HD models have a label beginning with H. The 

SECOND LETTER DENOTES THE INITIAL ROTATION PROFILE. THE LABELS OF BASIC RUNS WITH THE INITIAL RIGID ROTATION CONTAINS THE LETTER R. 

Models started with a solar-type or an anti-solar type DR have a letter S or A. The number indicates the given rotation rate 

El 0 WITH THE DECIMAL POINT REMOVED. THE MEAN QUANTITIES V rm s, #eq, Ro, T c , e M , AND e K ARE EVALUATED AT THE SATURATED STATE IN THE 
CONVECTIVE ENVELOPE. THE SUPER-CRITICALITY LEVEL S IS DEFINED BY (Ra - Ra c ) / Ra c , WHERE Ra c ~ 9.83 E" 116 W I TH THE EKMAN NUMBER OF 
E = v/(Qod$). Here we adopt, as the critical Rayleigh number, the fitting formula derived by|Al-Shamali et al. (2004j), in which 

THEY OBTAINED IT FROM THE SIMULATION OF THE SPHERICAL SHELL CONVECTION WITH VARYING THE SHELL THICKNESS AND THE EKMAN NUMBER. 
a e = V(f)(r = R,6 = 7t/2)/(Q 0 R) IS THE DR PARAMETER. In THE FINAL COLUMN, THE RESULTANT DR PROFILES, S (SOLAR-TYPE) OR AS (ANTI-SOLAR 

TYPE), ARE SUMMARIZED. 



Qq 

S 

Trms 

Beq 

T c 


£m 

Ro 

ROconv 

Q!e 

type 

MR00625 

0.0625 

155 

0.037 

0.028 

8.0 

1.3 x 10" 3 

2.7 x 10- 4 

6.272 

3.737 

-0.351 

AS 

MR0125 

0.125 

69 

0.040 

0.028 

7.6 

1.6 x 10" 3 

4.0 x 10- 4 

3.310 

1.868 

-0.386 

AS 

MR025 

0.25 

30 

0.029 

0.021 

10.2 

8.4 x KT 4 

8.8 x 10“ 4 

1.230 

0.934 

-0.094 

AS 

MR030 

0.30 

24 

0.028 

0.021 

10.5 

8.3 x lCT 4 

7.4 x 10“ 4 

0.995 

0.778 

-0.069 

AS 

MR035 

0.35 

20 

0.030 

0.022 

10.0 

8.6 x lCT 4 

2.6 x 10- 4 

0.898 

0.667 

0.045 

S 

MR040 

0.40 

17 

0.029 

0.021 

10.3 

8.3 x lCT 4 

2.4 x 10- 4 

0.763 

0.584 

0.086 

S 

MR050 

0.50 

13 

0.028 

0.020 

10.8 

8.0 x lO" 4 

2.5 x 10- 4 

0.580 

0.467 

0.088 

S 

MR060 

0.60 

10 

0.027 

0.020 

11.3 

7.4 x KT 4 

2.7 x 10“ 4 

0.463 

0.389 

0.078 

S 

MR150 

1.50 

3 

0.018 

0.012 

17.0 

2.8 x 1CT 4 

2.9 x 10“ 4 

0.123 

0.156 

0.011 

S 

MR300 

3.00 

0.75 

0.010 

0.007 

29.3 

9.1 x 1CT 5 

4.4 x 10- 4 

0.036 

0.078 

0.0003 

S 

HR00625 

0.0625 

155 

0.040 

0.029 

7.4 

1.5 x 1CT 3 

- 

6.761 

3.737 

-0.257 

AS 

HR025 

0.25 

30 

0.055 

0.040 

5.4 

4.0 x 10" 3 

- 

2.306 

0.934 

-0.384 

AS 

HR030 

0.30 

24 

0.056 

0.040 

5.4 

4.5 x 10“ 3 

- 

1.955 

0.778 

-0.371 

AS 

HR035 

0.35 

20 

0.056 

0.040 

5.3 

4.9 x 1CT 3 

- 

1.683 

0.667 

-0.310 

AS 

HR040 

0.40 

17 

0.057 

0.040 

5.3 

5.2 x 1CT 3 

- 

1.481 

0.584 

-0.296 

AS 

HR045 

0.45 

15 

0.034 

0.024 

00 

bo 

1.4 x 1CT 3 

- 

0.796 

0.519 

0.105 

S 

HR050 

0.50 

13 

0.034 

0.024 

8.9 

1.4 x 10" 3 

- 

0.707 

0.467 

0.114 

S 

HR055 

0.55 

12 

0.034 

0.024 

00 

bo 

1.5 x 10“ 3 

- 

0.649 

0.425 

0.135 

S 

HR060 

0.60 

10 

0.034 

0.024 

00 

bo 

1.6 x 10“ 3 

- 

0.593 

0.389 

0.132 

S 

HR300 

3.00 

0.75 

0.011 

0.009 

28.3 

2.2 x 1CT 4 

- 

0.037 

0.078 

0.010 

S 

MS025 

0.25 

30 

0.030 

0.021 

10.0 

8.4 x 1CT 4 

9.2 x 1CT 4 

1.254 

0.934 

-0.071 

AS 

MS030 

0.30 

24 

0.030 

0.021 

10.2 

7.7 x 1CT 4 

8.5 x 10- 4 

1.031 

0.778 

-0.050 

AS 

MS035 

0.35 

20 

0.030 

0.022 

10.0 

8.3 x 1CT 3 

2.9 x 10 ~ 4 

0.896 

0.667 

0.065 

S 

HS030 

0.30 

24 

0.053 

0.042 

5.7 

4.8 x 10" 3 

- 

1.837 

0.778 

-0.370 

AS 

HS035 

0.35 

20 

0.049 

0.042 

6.1 

5.0 x 10“ 3 

- 

1.462 

0.667 

-0.306 

AS 

HS040 

0.40 

17 

0.035 

0.026 

8.5 

1.4 x 10“ 3 

- 

0.928 

0.584 

0.102 

S 

MA030 

0.30 

24 

0.030 

0.021 

10.1 

7.7 x KT 4 

7.6 x 10“ 4 

1.034 

0.778 

-0.062 

AS 

MA035 

0.35 

20 

0.030 

0.022 

9.9 

9.3 x 1CT 4 

4.0 x 1CT 4 

0.909 

0.667 

-0.017 

AS 

MA040 

0.40 

17 

0.030 

0.021 

10.1 

7.8 x 1CT 4 

2.7 x 1CT 4 

0.780 

0.584 

0.059 

S 

HA040 

0.40 

17 

0.056 

0.040 

5.4 

5.3 x lO" 3 

- 

1.459 

0.584 

-0.293 

AS 

HA045 

0.45 

15 

0.055 

0.040 

5.4 

5.3 x 10" 3 

- 

1.286 

0.519 

-0.281 

AS 

HA050 

0.50 

13 

0.046 

0.024 

6.5 

1.4 x 10“ 3 

- 

0.974 

0.467 

0.111 

S 


convective eddies (this can be justified because of our weakly 
stratified system), = [(3/2)(i^ + vj ;)] 1 / 2 is the mean ve¬ 
locity, and Ta = (2f] 0 d 2 /z/) 2 is the Taylor number. Angular 
brackets denote time- and volume-average “in the convective 
envelope" at the saturated state. The definitions of v rrns and 
Co are almost identical to those in Kapyl a et al.| ( |2Ql l[ [2014) 
and |Karak et al.| ( |2014| ). Note that, in this study, the convec- 
tive Rossby numbeFis only a function of the rotation rate, i.e., 
Roconv Stf, because the other parameters are fixed in our 
simulations (see following paragraphs). 

The stress-free boundary condition for the velocity is im¬ 
posed on both the radial boundaries. As for the magnetic field, 
we assume the perfect conductor at the inner boundary and the 
radial field condition at the outer boundary. A constant energy 
flux which drives the convective motion is imposed at the in¬ 
ner boundary and the temperature is fixed to be T s at the outer 
boundary. 

The equations Q -0 are solved by the second-order cen¬ 
tral difference s cheme with spatial discr etization with Yin- 
Yang grid. See |Kageyama & Sato| ( |2004| and MYK13 for 
details. Non-dimensional parameters Pr = 0.27, Pm = 4.0, 


and Ra = 2.2 x 10 5 (i.e., p = 7.8x 1CT 5 , rj = 2.0 x 10 -5 , and 
k, = 3.9 x 1CT 4 ) are commonly adopted in all the models stud¬ 
ied here. The total grid size for all the simulation runs is 81 
(in r) x 131 (in 6) x393 (in <fi) x2 (Yin & Yang). To explore 
the response of the large-scale flow to the rotation rate, we 
vary the magnitude of Qq in the MHD and HD models while 
keeping the background hydrostatic state unchanged. A ran¬ 
dom temperature perturbation and small “seed” magnetic field 
(only for MHD models) are introduced at the same time in the 
convection zone when the calculation starts. 

3. NUMERICAL RESULTS 

Our simulation models are summarized in Table 1. Model 
names are given in the left most column. The first letter in the 
model name denotes the presence and absence of the magnetic 
field. The labels of MHD models contain the letter M. HD 
models have a label beginning with H. The second letter de¬ 
notes the initial rotation profile. The labels of basic runs with 
the initial rigid rotation contains the letter R. Models started 
with a solar-type or an anti-solar type DR have a letter S or 
A. The number indicates the given rotation rate Dq with the 
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Fig. 1.— The radial velocity profile v r (r, 6) at r = 0.85 R when t ~ 360r c on a spherical surface. The normalization unit is the mean velocity v rm s of each model. 
The panels (a)-(c) correspond to the MHD models with Qq = 0-3, 0.35 and 0.4, and the panels (d)-(f) are for the HD models with Qq = 0-4, 0.45 and 0.5. The 
lighter and darker tones depict upflow and downflow velocities. 


decimal point removed. 

The mean quantities evaluated at the saturated state (v rms , 
B qq , t c , ck, and ck), the Rossby number (Ro), and the 
convective Rossby number (Ro conv ) are shown in Table 1, 
where r c = d /v rms is the convective turn-over time, B Qq = 
((3/2 )p (v^ + v 2 )) 1 / 2 is the the equipartition strength of mag¬ 
netic field, ck = (pf 2 )f s /2 and cm = {B 2 ) f s /2 are the mean 
kinetic and magnetic energy densities. Angular brackets with 
subscript “fs" denote time- and volume-average in the full- 
spherical shell domain at the saturated state. As a measure of 
the convection strength, the estimated super-criticality level, 
defined by S = (Ra-Ra c )/Ra c , is shown in the second column, 
where Ra c is the critical Rayleigh number (see caption for the 
definition of Ra c ). In the second last and last columns, DR 
parameter a e (see §3.2 for the definition) and, the resultant 
DR profile, S (solar-type) or AS (anti-solar type), are summa¬ 
rized. 

We perform three sets of runs. In the first set, we run models 
from the rigid rotation (Set R: Runs MR and HR) with the 
initial conditions described in §2. In the second and third sets, 
we examine the bistability of the rotation profile by taking 
either a solar-type (Set S: Runs MS and HS) and an anti-solar 
type (Set A: Runs MA and HA) solution as initial conditions. 
The MHD effects on the convective flows are studied in each 
set. 


3.1. Convective Flows in Basic Runs with Initial Rigid 
Rotation 

After the convective motion sets in, it reaches a saturated 
state before ~ 100 r c in the basic run with the initial rigid ro¬ 
tation. We have run the simulations up to around 400r c , which 
is almost comparable to the magnetic diffusion time, and stud¬ 
ied the properties of the convective flow at the equilibrated 
state. 

Figure 1 shows, in the Mollweide projection, the radial 
velocity distribution when t ~ 360r c on spherical surface at 
r = 0.85R. The normalization unit is the mean velocity, Vnns, 
of each model. Panels (a)-(c) correspond to the MHD mod¬ 
els (MR030, MR040, and MR045) and panels (d)-(f) are for 
the HD models (HR040, HR045, and HR050). The lighter 
and darker tones depict upflow and downflow velocities in the 
range |v r /vn„ s | < 0.5. 


The rotating stratified convection is characterized by up- 
down asymmetry in all the models: the slower upflow dom¬ 
inant cells surrounded by the faster and narro wer downflow 
lanes (e.g., |Spruit et aL 1990[ |Miesch||2005| ). As the rota¬ 
tion rate increases, the convective cell shrinks and its structure 
is changed from the cellular pattern to the elongated colum 
nar pattern formed in the equatorial region (e.g. jBusse 1970 

■ * ” *or ” 


Glatzmaier & Gilman 1981; Brummell et al. 1996, 1998; Mi- 


esch et al. 2000]). The transition from the cellular to c olumnar 

convections has been obsered in earlier simulation s (Kapyla 


|2011||20T4l|Gastine et al. 2013, 2014] |Guerrero et aT 

2013) and occurs at around Do = 0.35 (Ro con v = 0.67) in the 


et al. 


MHD model but D 0 = 0.45 (Ro conv = 0.52) in the HD model. 
The magnetic field affects the convective motion, suggesting 
its impacts on the turbulent angular momentum transport and 
thus the resultant large-scale DR. 

Shown in Figure 2 is the profile of the mean DR, defined 
by D(r,$) = V 0 /(rsin#) + Do, where the overbar denotes the 
time- and azimuthal-average. The normalization unit is D 0 
of each model. The time average spans in the range of 
200 < t/r c < 300. Panels (a)-(c) correspond to the MHD 
models and panels (d)-(f) are the HD models. The red (blue) 
tone denotes the higher (lower) angular velocity. The white 
lines are iso-rotation contours. The interface between the con¬ 
vective envelope and the radiative layer is denoted by black 
dashed line. 

As shown in earlier studies, there exists two regimes in the 
rotation profile: an anti-solar type DR with slower equator 
[(a) & (d)] and a solar-type DR with faster equator [(b),(c),(e) 
& (f)j (e.g., IGilmanl 19771 |K apyla et al.|2011 1 |20l?l |Gastine 

120141 IGuerrero et 


|et al |2Q \3\ |20 1 4[ |Guerrero et al.|2013| ). While the anti-solar 

type DR is established in the model with the smaller rotation 
rate, the solar-type DR develops in the regime of the higher 
rotation rate. The most remarkable difference between the 
MHD and HD models is the critical rotation rate (and thus the 
critical convective Rossby number) for the transition. While 
the transition between the anti-solar and solar-type DR pro¬ 
files occurs at around D 0 = 0.35 (Ro conv = 0.67) in the MHD 
model, it does at around Do = 0.45 (Ro con v = 0.52) in the HD 
model. 

The DR profile is dominated by so-called “Taylor- 
Proudman balance" both in the MHD and HD models (e.g., 
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Fig. 2.— Time-averaged mean angular velocity 6). The panels (a)-(c) correspond to the MHD models with Qq = 0-3, 0.35 and 0.4, and the panels (d)-(f) 
are for the HD models with Qq = 0.4, 0.45 and 0.5. The white solid line denotes the iso-rotation contour. The black dashed line denotes the interface between 
convective and radiative zones. 
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|Pedlosky |1987| ) and is still far from the solar conical isotachs 
deduced from the helioseismic measure ment (e.g., |Kitchati 
nov & Rudiger |1995[ |Miesch||2005 |Rempel|j2QQ5 
| |2011| ). Nevertheless, it is intriguing that the is( 


Masada 


iso-rotation 

contour is more conical in the MHD model than in the HD 
model. Especially, the radial shape of the isotachs is pro¬ 
nounced around the equator in the MHD model near the tran¬ 
sition (see MR035). 

When comparing the two models at the same rotation rate, 
we can see that the rotational shear is generally weaker in 
the MHD model. A strong polar jet (polar vortex) appeared 
in the HD model is also suppressed in the MHD model (see 
HR045 and HR050). In addition, the MHD effect seems to 
suppress the (viscous) spreading of the rotation profile of the 
convective envelope into the radiative layer in all the MHD 
models, as was observed in §4 of MY K13 (see, e.g., Spiegel 
& Zahn|l992[ [Gough & Mclntyre|1998| ). The reduction of the 


rotational shear by the influence of the magnetic field s agrees 
with that has been reported in the earlier studies (e.g.,|Brun et 


al.|2004[ [Beaudoin et al.|2013] Fan & Fang|2014||Karak et ah 


201 4|) and is discussed further in §3.2. 

Figure 3 shows the time-averaged mean meridional flows. 
Panels (a)-(c) correspond to the MHD models and panels (d)- 
(f) are the HD models around the transition between the solar 
and anti-solar type rotation profiles. The color contour depicts 
the meridional flow velocity, defined by v m = (vj; + v^) 1 / 2 . The 
normalization unit is of each model. The streamlines are 
over-plotted with a length proportional to the flow speed. The 
white dashed line denotes the interface between the convec¬ 
tive and radiative layers. 

The circulation flow is primarily counter-clockwise (clock¬ 
wise) in the bulk of the convection zone in the northern (south¬ 
ern) hemisphere. This is consistent with the observati ons of 
the solar sub-surface meridional flow (e.g.JKomm et al.[ 1993 


|Giles et al.|1997 |UMch|2010||Zhao et al.|2013| ), and a com 

mon feature of all the models. However, there is a difference 
in the circulation pattern. While a large single-cell is formed 
in the model with the anti-solar type DR (see MR030 and 
HR040), the model with the solar-type DR shows a multiple¬ 
cell pattern (see MR040 and HR050). A similar transition 
from single to multiple cell patterns of the meridional flow has 
been observed before in different settings (e.g.,|Kapyla et al. 


|2011[|2014HGastine et al.|2013[|Karak et al.|2014| ). We note 

that, as well as the rotation pattern, the transition of the merid¬ 
ional flow pattern occurs at lower (thus higher Ro conv ) in 
the MHD model than in the HD model. 

3.2. Control Parameter of Differential Rotation 

The transition between the anti-solar and solar-type DRs 
is a result of the structural difference of the whole three- 
dimensional convective motion separated at the critical rota¬ 
tion rate. Here we quantitatively examine the dependences of 
the DR profiles on two key diagnostic parameters, the convec¬ 
tive Rossby number Ro conv and the Rossby number Ro. Then 
we address which can better capture the transition between 
the anti-solar type and solar-type DRs. Following Gastine et 
al. ( 2014| , we quantify the DR by the amplitude of the mean 
azimuthal flow at the equatorial surface: 


OLp 


v < f,(r = R,0 = 7t/2) 


q 0 r 


( 8 ) 


The positive and negative a e denote the solar-type and anti¬ 
solar type DRs. 


Figure 4(a) shows the DR parameter (a e ) as functions of the 
convective Rossby number (Ro con v) [bottom axis] and the cor¬ 
responding rotation rate (Qq) [top axis]. The blue crosses and 
red diamonds denote the MHD and HD models. The horizon¬ 
tal and vertical dashed lines indicate a e = 0 and Ro conv = 1. 
As stated in § 3.1, the DR profile transits from the solar to 
anti-solar type regimes as the convective Rossby number in¬ 
creases. The absolute amplitude of a e is, on average, larger 
in the anti-solar type DR than in the solar-type DR, and takes 
the maximum and minimum around the transition. These are 
common featu res of the two models and are compatible with 
earlier stu dies (Ga stine et al.|2014[|Kapyla et al.|2014[|Karak 
|et al.||2014| ). The critical value of Ro conv for the transition is 
different in the two models. It occurs when 0.67 < Ro con v < 
0.78 for the MHD model and 0.52 < Ro conv < 0-58 for the 
HD model, indicating a more gradual transition under the in¬ 
fluence of the magnetic field. When comparing two models at 
the same Ro conv > the magnitude of a e is smaller in the MHD 
model than in the HD model. We confirm that the magnetic 
field has a role in suppressing the DR and helps to produce 
solar-type DR, whi ch have be en suggested in MYK13, [Fan &| 
Fang (2014| andjKarak et al. (|2014|). 


Uur intriguing finding is that the critical value of the Rossby 
number for the transition is almost the same between the 
MHD and HD models. Shown in Figure 4(b) is the depen¬ 
dence of the DR parameter on the Rossby number (Ro). The 
blue crosses and red diamonds are corresponding to the MHD 
and HD models. The horizontal and vertical dashed lines de¬ 
note a e = 0 and Ro = 1. The DR parameter decreases with 
the increase of the Rossby number in both models around the 
transition. The transition between the solar and anti-solar type 
DR profiles occurs at Ro ~ 1.0 both in the MHD and HD mod¬ 
els. Not only the critical value, the sharpness of the transition 
is also similar between the two models. The critical Rossby 
number for the transition obtained here gives a good agree- 
ment with that obta i ned in the hydrodynamic simulations by 
|Kapyla et al.| ( [201 1| |2014| ). These indicate that the transition 
between the two rotation regimes is controlled by the Rossby 
number rather than by the convective Rossby number (i.e., the 
rotation rate itself in our study) regardless of the presence of 
the magnetic field. 

We stress here thatfKarak et al.|(|2014|) also studied the mag¬ 
netic influence on the DR profile by comparing HD and MHD 
models, but did not discuss their Ro-dependences. Reading 
the critical Rossby num ber for the tra nsition from Table 1 of 
Karak et al. ( |2Q14| ) and |Kapyla et al.| ( |2014| , it is found that 


the transition occurs at a similar value of Ro in their MHD and 
HD models as well . Our finding is thus consistent with the 
earlier studies. The influence of the magnetic field on the DR 
profiles was also briefly mentioned in Gastine et al. ( 2014| ). 
However, they focused on the dependence of the DR profile 
on the “local Rossby number (Ro/)" and did not examine its 
Ro-dependence. Here the local Rossby number is defined by 
R 0 / = l u 5u Tms /(jtQ,Qd), where Su Tms is the RMS value of the 
non-axisymmetric (fluctuating) component of the velocity (in¬ 
cluding both meridional and azimuthal components) and l u is 
the mean spherical harmonic degree obtained from the kinetic 
energy spectrum (e.g., [Christensen & Aubert||||2006[ |Schrin 


ner et al.||2012| ). Since they showed that the critical “local' 


Rossby number for the transition is different between MHD 
and HD models, we can say that the Rossby number is even 
better than the local Rossby number at least to describe the 
transition between the anti-solar and solar-type DR profiles in 
a unified manner. 
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Fig. 3.— Time-averaged mean meridional flows. Panels (a)-(c) correspond to the MHD models (MR030, MR035, MR040) and panels (d)-(f) are the HD 
models (HR040, HR045, HR050). The color contour depicts the meridional flow velocity, defined by v m = [vjr + v^] 1 / 2 normalized by Vrms of each model. The 
streamlines are overplotted with a length proportional to the flow speed. The white dashed line denotes the interface between convection and radiative zones. 
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In order to reveal a cause of the difference in the critical 
value of the convective Rossby number between the two mod¬ 
els, we present, in Figure 5, the dependence of the mean con¬ 
vective velocity (v rms ) on the convective Rossby number for 
the MHD (blue crosses) and HD (red diamonds) models. Here 
the regime around the transition is focused. When compar¬ 
ing two models at the same Ro con v, the MHD model has the 
smaller convective velocity than the HD model. The MHD 
model thus can maintain the solar-type DR even at the lower 
rotation rate (thus higher convective Rossby number) in com¬ 
parison with the HD model. Since the Lorentz force of the 
dynamo-generated magnetic field puts the brakes on the con - 
vective motion (e.g., |Cattaneo et al.|2003[|Brun et al.|2004| , 
the MHD model becomes more rotationally-constrained than 
the HD model. This could be the reason why the magnetic 
field helps to produce the solar-type DR. |Fan & Fang ( 2014) 
and |Karak et ah] ( |2014| ) also argued that the reduction of the 
convective velocity Fy the dynamo-generated magnetic field 
is responsible for the difference of the critical rotation rate 
(i.e., critical convective Rossby number) for the solar to anti¬ 
solar transition between the HD and MHD models. It is in¬ 
teresting that such a drastic difference between the MHD and 
HD models observed in (Fig.5) is unified in the view of 
the Ro-dependence (Fig.4b). 

3.3. Bistability of the Rotation Profile 

The transition between the solar and anti-solar type rota¬ 
tion patterns has been known since |Gilman|(|1977|) and stud¬ 
ied in various setting s ([Kapyla et al.|201 if Gastine et al.|2013[ 
Guer rero et al.|2013] ). However, the bistability of the rotation 
profile has not been sufficiently -studied because it has only 
recen tly been discovered ( Gasti ne et al.||2014} |Kapyla et al. 
|2014| ). We revisit here the dependence of the rotation profile 
on the history both in the MHD and HD models by taking ei¬ 
ther a solar-type (Set S: Runs MS and HS) and an anti-solar 
type (Set A: Runs MA and HA) solution as initial conditions. 

In the Set S, we perform simulations by starting from the 
saturated state of MR040 (HR050) with the solar-type DR for 
the MHD (HD) run and decrease the rotation rate. In contrast, 
in the Set A, we choose the final state of MR025 (HR030) 
with the anti-solar type DR as the progenitor for the MHD 
(HD) model and increase the rotation rate. The DR parameter 
(a e ) for these models are shown in Figure 6(a) as a function 
of the convective Rossby number (Ro con v)- The models MS, 
MA, HS and HA are denoted by cyan triangles, green circles, 
pink crosses and orange squares, respectively. The DR pa¬ 
rameters of the basic runs (MR and HR) are also shown by 
blue crosses and red diamond s as referen ces. 

As shown in |Gastine et al.| ( [20T4] ) and |Kapyla et al.| ( |2014| ), 
the bistability of the rotation profile can be observed near the 
transition in the HD model. In the Set HS, the transition of 
the HD solution occurs when 0.58 < Ro con v <0.67 while it 
does when 0.47 < Ro con v < 0-52 in the Set HA. On the other 
hand, the set MS started from the solar-type DR reproduces 
the similar results with the set MR with the initial rigid rota¬ 
tion (see cyan triangles an d blue crosse s ). Thi s is consistent 
wit h the r ecent studies of |Fan & Fang| ( [20l4] ) and |Karak et 
[ah] ([20 14) who found in their convective dynamo simulations 
mat the bistability is disappeared and the DR profile becomes 
independent from the hysteresis when allowing the growth of 
the MHD dynamo. However, in the set MA which is started 
from the anti-solar type DR, the bistability seems to be not 
disappeared (see green circles). 

The remaining of the bistability in the set MA of our MHD 


models may be related to the presence of the stable layer un¬ 
derlying the convective envelope, which is a difference be¬ 
tween our simulati on model and the m odels adopted in |Fan & 
Fang | (|2014|) andlKarak et al.|(|2014|). Carefully studying the 
rotation profile of the MR03T) (Fig.2a) which has a very si 


profile of the MR030 (Fig.2a) which has a very simi¬ 
lar rotation profile with the progenitor for the set MA, we can 
find that the stable layer has a strong equatorial acceleration 
despite the anti-solar type DR established in the convection 
zone. Since the differential rotation in the stable layer should 
change the latitudinal temperature and entropy distributions 
to achieve the thermal wind ba lance there (e.g., |Kitchatinov| 
& Rudige r| 1995 [ [Rempel 2005), it may enforce the anti-solar 
type DR even in the high i 2o regime which should provide 
the solar-type DRs. Our speculation on the role of the sta¬ 
ble layer in the maintenance of the internal rotation will be 
quantitatively examined in a subsequent paper. 

It is interesting that, by using the Rossby number (Ro), we 
can again unify the transition between the solar and anti-solar 
type rotation profiles independently of the hysteresis. Shown 
in Figure 6(b) is the DR parameters of all the models as a 
function of the Rossby number. The symbols have the same 
meanings as those in Figure 6(a). It is found that the tran¬ 
sition occurs at Ro ~ 1 in all the simulation sets. Not only 
the critical value, the sharpness of the transition is also sim¬ 
ilar between all the cases. This confirms that the transition 
between the two rotation regimes is better captured by the 
Rossby number rather than by the convective Rossby number. 

One important finding in our study is that the initial rotation 
profile, i.e., the evolution history of the stellar rotation, has an 
impact on the convective velocity if we does not allow the 
growth of the magnetic field. The dependence of the mean 
convective velocity (v^s) on the convective Rossby number 
(Ro conv ) is shown in Figure 7. The symbols have the same 
meanings as those in Figure 6(a). Here the regime around the 
transition is focused. The bistable nature of the convective ve¬ 
locity can be seen in the HD models, while the MHD models 
have similar convective velocity around the transition. Thus, 
in the HD cases, the Set S (Set A) with the initial solar-type 
(anti-solar type) DR becomes more rotationally-dominated 
(inertia dominated) than the other sets. In contrast, in the 
MHD cases, the evolution history of the DR does not change 
the convective velocity and the rotational dominance in the 
system. This would be due to the regulation of the convec¬ 
tive velocity by the dynamo-generated magnetic field. These 
indicate that the bistability appeared in the rotation profile is 
essentially a consequence of the hysteresis of the convective 
velocity, which can be weakened by the dynamo-generated 
magnetic field. 

Since the convective Rossby number denotes the relative 
importance of the buoyancy force in comparison with the 
Coriolis force, the magnetic effect and the evolution history 
are not directly reflected in its definition. In contrast, the 
Rossby number describes the significance of the inertia force 
arising as the results of the hysteresis and the non-linear inter¬ 
action between the convective motion and the magnetic field, 
relative to the Coriolis force. Therefore, the effects of the 
magnetic field and the hysteresis are naturally reflected in it. 
This would be the reason why the transition between the solar 
and anti-solar type rotation profiles can be unified in the view 
of the Ro-dependence. 


3.4. Dependence of Magnetic Dynamo Activity on Rotation 

Rate 
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Fig. 4.-— (a) DR parameter a e as functions of the convective Rossby number (Ro con v) and the corresponding rotation rate (Qq)- The horizontal and vertical 
dashed lines indicate a e = 0 and Ro CO nv = 1. (b) DR parameter a e as a function of the Rossby number (Ro). The horizontal and vertical dashed lines indicate 
a e = 0 and Ro = 1. The blue crosses and red diamonds denote the MHD and HD models with the initial rigid rotation in both panels. 



Alu conv 

Fig. 5.— Mean convective velocity v rm s as a function of the convective 
Rossby number (Ro CO nv)- The blue crosses and red diamonds denote the 
MHD and HD models with the initial rigid rotation. 

The MHD model exhibits a rich variety of the magnetic 
dynamo activity depending on the rotation rate. We show 
the time-latitude diagram of (a )(B r ) < f > (left column), (b )(Be) < f ) 
(middle column) and (c)(#</>)</> (right column) of MR00625, 
025, 030, 050, 060, and 150 at r - 0.65R (mid-stable zone) 
in Figure 8 and at r- 0.85R (mid-convection zone) in Figure 
9. Here the angular brackets with subscript <fi denote the az¬ 
imuthal average at a given depth. The time and amplitude are 
normalized by r c and B eq , respectively. The red and blue tones 
denote the positive and negative mean-field strength. The hor¬ 


izontal dashed-line in each panel denotes the equator. The 
upper three models have the anti-solar type DR and the lower 
three models are characterized by the solar-type DR. MR025- 
060 are the models near the transition. 

The model with the smallest rotation rate is dominated 
by the turbulent magnetic field and has little large-scale 
component both in the radiative and convective layers (see 
MR00625). This would be due to the weak Coriolis force, 
yielding weak differential rotation. At around the transition of 
the rotation profile, the spatiotemporal coherence of the mag¬ 
netic field becomes remarkable. In the slow rotation regime 
with the anti-solar type DR, a quasi-steady large-scale mag¬ 
netic field with a dipole symmetry is organized (see MR025 
and 030). (5 r )^ and (Bq)^ are stronger in the convection 
zone rather than in the radiative zone, while (#</>)</> is stronger 
in the stable zone than in the convection zone. Among the 
three magnetic components, (B^)^ is dominant in these mod¬ 
els. The large-scale B ^-component is the strongest at the mid¬ 
latitude and is antisymmetric about the equator. This is con¬ 
sistent with the shearing of the Bq -component, i.e., (^-process, 
by the anti-solar type DR. 

The properties of the magnetic dynamo change drastically 
across the transition of the DR profile. The large-scale mag¬ 
netic field shows the polarity reversal in the fast rotation 
regime with the solar-type DR. In the models near the tran¬ 
sition, i.e., MR050 and 060, the large-scale magnetic com¬ 
ponent with polarity reversals is built up mainly in the stable 
zone. It has commonly B^ -dominance. The large-scale B$- 
component is the strongest at around the equator in these mod¬ 
els unlike the models with the anti-solar type DR in which it 
is the strongest at the mid-latitude. In MR050, the convection 
zone is dominated by the turbulent magnetic field. However, 
we can observe a weak spatiotemporal coherence of the mag¬ 
netic field in the convection zone in MR060. It is interesting 
that the cycle of the polarity reversal of the large-scale mag- 
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Fig. 6. — DR parameters ( a e ) of the models with non-rigid initial rotations as functions of (a) the convective Rossby number (Ro CO nv) and (b) the Rossby 
number (Ro). The models MS, MA, HS and HA are denoted by cyan triangles, green circles, pink crosses and orange squares, respectively in both panels. The 
horizontal and vertical dashed lines indicate (a) a e = 0 & Ro CO nv = 1 and (b) a e = 0 and Ro = 1. The DR parameters of the basic runs (MR and HR) are shown by 
blue crosses and red diamonds as references. 



conv 

Fig. 7.— Mean convective velocity v rm s as a function of the convective 
Rossby number (Ro con v)- The models MR, MS, MA, HR, HS and HA are 
denoted by blue crosses, cyan triangles, green circles, red diamonds, pink 
crosses and orange squares, respectively. 


netic field around the equator is shorter in the convection zone 
than in the radiative zone. 

The most intriguing dynamo activity can be observed in the 
model with a fast rotation beyond the transition (see MR 150). 
Note that this model has 3-4 times higher rotation rate than 
that adopted in the model around the transition. In the sta¬ 
ble zone, the strong large-scale magnetic component with a 
dipole symmetry is built up. It concentrates at the high latitu¬ 
dinal region (~ 80°) and seems to have a long cycle period of 
the polarity reversal (r cyc > 200r c ). The spatial distribution of 


the -field suggests that the fTeffect plays a crucial role in 
amplifying the large-scale magnetic field in the stable region. 
In the convective envelope, two types of the dynamo mode can 
be clearly observed. The long-lasting magnetic component 
with a dipole symmetry, which is similar to that in the stable 
zone, is found in the high latitude. In contrast, at around the 
equatorial region, the oscillatory dynamo mode with a shorter 
cycle period (r cyc ~ 10r c ) and a poleward migration is excited 
and sustained. The difference of the dynamo mechanism be¬ 
tween the models with the different rotation profiles is dis¬ 
cussed in §4. 

4. DISCUSSION: THE RELATION BETWEEN MAGNETIC DYNAMO 
ACTIVITY AND KINETIC HELICITY 

The purpose of this paper is to elucidate the effect of the 
magnetic field on the formation of the rotation profile. Nev¬ 
ertheless, a rich variety of the magnetic dynamo activity ob¬ 
served in the MHD models (see §3.4) makes us interested in 
their dynamo mechanisms. Here we discuss a possible key 
ingredient for the oscillatory magnetic dynamo with focusing 
on its relation to the kinetic helicity of the convection flow. 

As an indicator of the magnetic dynamo activity, we study 
the relative amplitude of the magnetic energy to the kinetic 
energy which is measured by the ratio cm/ck, where cm and 
ck are mean magnetic and kinetic energy densities at the sat¬ 
urated state (see §3 for the definitions of cm and ck). In ad¬ 
dition, the mean kinetic helicity (net helicity), H av e> of each 
MHD model is also evaluated, where it is defined by 

n^ = (\n N \+\n s \)/2, (9) 

with 

Hi s (v • (V x v))i (/ = N,S). (10) 

Here the angular brackets with subscripts “N" and “S" denote 
the time- and volume-average in the northern and southern 
hemispheric convection zones at the saturated state. 
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Fig. 8. — Time-latitude diagram of (a) ( B r (left column), (b) ( Bq (middle column) and (c) (#</>)</> (right column) at r = 0.65 R (mid-stable zone) for the 
models MR00625, 025, 030, 050, 060, and 150, respectively. The time and the amplitude are normalized by r c and 5 eq . The red and blue tones denote the 
positive and negative mean-field strength. Horizontal black dashed-line corresponds to the equator. The top three models have the anti-solar type DR and the 
bottom three models have the solar-type DR. 


Shown in Figure 10 is cm/ck (blue line with diamonds 
[right axis]) and H av e (fed line with circles [left axis]) as a 
function of the Rossby number (Ro) for the basic MHD runs 
(Set MR). The vertical dashed line denotes Ro = 1. The ratio 
cm/£k has a bit complicated profile: In the regime of the small 
or large Rossby number, that is Ro <C 1 or Ro 1, cm/^k de¬ 
creases with the increase of Ro. In contrast, at around the 
transition, the higher Ro rather provides the larger 6 m/ck- 
The mean kinetic helicity shows a single peak profile with 
the maximum at around Ro ~ 1 where the transition from the 
solar to anti-solar type rotation profiles occurs. In the regime 
of Ro < 1 where the solar-type DR is established, €m/ck anti¬ 
correlates with Have- By contrast, the ratio decreases with the 
kinetic helicity in the anti-solar type rotation regime of Ro > 
1. Since the temporal property of the large-scale magnetic 


field also changes from the oscillatory to the stationary across 
Ro ~ 1 as was described in §3.4, this suggests that the dynamo 
mechanism would be associated in some way with the kinetic 
helicity in the convection zone. 

The Ro-dependence of the net kinetic helicity, which shows 
a peak at moderate rotation rate, is deeply related to the for¬ 
mation of the convective column. Since the angular velocity 
is lower in the regime of the anti-solar type DR or Ro ^ 1, 
the up-down asymmetry of the convective motion is lower, 
leading to less preference for one sign of the kinetic helic¬ 
ity. On the other hand, in the regime of the solar-type DR or 
Ro <C 1, the convection has a tendency to align in columnar 
rolls parallel to the rotation axis (see §3.1). In the convection 
columns, the velocity is mostly perpendicular to the rotation 
axis while the vorticity is mostly parallel to the rotation axis, 
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Fig. 9. — Same as Figure 8 but at r = 0.85 R (mid-convection zone). 


resulting in a small net kinetic helicity (e.g., 


1981[ |Bro wning 2008). It was also clearly demonstrated by 


Knobloch et al. 


Sprague et al. (2006) with asymptotic numerical analysis of 
unstratified turbulence that the largest net kinetic helicity is 
established at moderate rotation rate and decreases as rotation 
becomes even more rapid (see also, |Kapyla et al.[2009| ). 

The anti-correlation between Ro and 6 M /ck in the solar- 
type DR regime would be consistent with the earlier studies 
(e.g., |Brown et al.||2011| |Kapyla et al.||2013| ). In addition, 
the correlation between Ro and 6m/£k around t he transition 


is also compatible with the numerical result of [Karak et al. 
p014| (see their Table 2). However, the role of the kinetic he¬ 
licity in the magnetic dynamo mechanism and its hidden con¬ 
nection to the observed relation ship between stellar rotation 


and magnetic activity (e.g., see |Gudel||2007t [Reiners 2012 
for reviews) remain to be fully elucidated although there are a 
great deal of work exploring the linkage between them (e.g., i>e ~ (ft) r = (v(Vxv)) r , is also shown in Figure 12, where 


Knobloch et al.|198l|pe nnings & We iss|199H Baliunas ^et al.| 

1996[ | Brandenburg et al.|1998||2012 ). 

To get a better grasp of the relation of the kinetic helic¬ 
ity and the magnetic dynamo mechanism, in Figure 11, we 
show the distribution of the time and azimuthally-averaged ki¬ 
netic helicity, %(r, 6) = v • (V x v), for the models MR00625, 
025, 030, 050, 060, and 150, respectively. The red and blue 
tones denote the positive and negative helicity. The horizon¬ 
tal black dotted-line corresponds to the equator. The dashed 
curve denotes the interface between the convection and ra¬ 
diative zones. The models listed in the upper three panels 
have the anti-solar type DR (MR00625, 025 and 030) and the 
lower three panels correspond to the models with the solar- 
type DR (MR050, 060 and 150). To clearly demonstrate the 
latitudinal distribution of the kinetic helicity, the radial av¬ 
erage of the time- and azimuthally-averaged kinetic helicity, 
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Fig. 10.— Mean helicity defined by "Have = (|7 ^n| + l^sD/2 (red circles) 
and ratio cm /£k (blue diamonds) as a function of the Rossby number (Ro) 
for the basic MHD runs. 


angular brackets with subscript r denote the radial average. 
The different line-types denote the models with different Do¬ 
lt would be trivial that the kinetic helicity has an antisym¬ 
metric profile with respect to the equator in the model with 
a sufficiently fast rotation. In the northern (southern) hemi¬ 
sphere, the negative (positive) helicity becomes predominant. 
This is because the downflow with the negative (positive) he¬ 
licity is on average stronger than the upflow with the posi¬ 
tive (negative) helicity in the northern (southern) hemisphere 
in the density stratified system (e.g., Spruit et ak] |1990| [Mi- 


esch||2005] ). The weaker equatorial antisymmetry of the kT 


netic helicity profile in the model with the smallest rotation 
(MR00625) would be a result of a weaker Coriolis force 
which provides less impact on the convective motion. The 
peak amplitude of the kinetic helicity is little dependent on 
the rotation rate (within factor 2). However, we can find from 
Figure 12 that the latitude for the peak of the kinetic helicity 
decreases as the rotation rate increases: While the higher lati¬ 
tudinal region has a larger helicity in the models with the anti¬ 
solar type DR (MR00625, 025, and 030), the lower latitudinal 
region has a larger helicity in the models with the solar-type 
DR (MR050, 060, and 150). 

It is the most remarkable that both the hemispheric and 
inter-hemispheric latitudinal gradients of the kinetic helicity 
become larger around the equator with the increase of the 
rotation rate. In the models with the solar-type DR and the 
oscillatory large-scale magnetic field, the kinetic helicity con¬ 
centrates on around the equatorial plane and has strong hemi¬ 
spheric and inter-hemispheric latitudinal gradients there. In 
contrast, the models with the anti-solar type DR and the sta¬ 
tionary large-scale magnetic field do not show the significant 
latitudinal gradient of the kinetic helicity around the equato¬ 
rial region. Since, when comparing Figures 9 and 11, the lat¬ 
itude providing the stronger large-scale magnetic field over¬ 
laps with that with the larger kinetic helicity, it is expected 
that the latitudinal distribution of the kinetic helicity is a key 
to distinguish the dynamo mechanisms. The larger the lat¬ 
itudinal gradient of the kinetic helicity around the equator 
becomes, the more the oscillatory large-scale magnetic field 
would be developed. 

Recently, [Mitra et al. ( |2010| ) reported an intriguing find¬ 


ing in their dynamo simulation by a forced helical turbulence 
in a wedge-shaped spherical shell that the oscillatory large- 
scale magnetic field with equatorial migration can be orga¬ 
nized, even without the D-effect, essentially due to an inter- 
hemispheric gradient of the kinetic helicity. They argued that 
the a 2 -dynamo mode excited by the helical turbulence with 
the equatorial antisymmetry is responsible for the oscilla¬ 
tory property. Furthermore, a connection between a 2 -dynamo 
mode and solar magnetism was discussed in some recentre- 
sults of convective MHD dynamos simulations (e.g.,|Simard| 
[et al.|2013{ [Masada & Sano|2014alb| ). 

When considering the latitudinal gradient of the kinetic he¬ 
licity around the equator as an important ingredient of the os¬ 
cillatory property of the dynamo even in our simulations, we 
can speculate the reason why there exists a difference of the 
dynamo property depending on the regime of the rotation pro¬ 
file. Since the model with the smaller rotation rate and thus 
the anti-solar type DR does not have a sufficient latitudinal 
gradient of the kinetic helicity around the equator to excite 
the oscillatory a 2 -dynamo mode, the D-effect solely works to 
build up the large-scale magnetic field and thus yields a sta¬ 
tionary solution. In contrast, the model with the faster rotation 
and thus the solar-type DR has a sufficient concentration and 
strong latitudinal gradient of the kinetic helicity around the 
equator, yielding the oscillatory large-scale magnetic compo¬ 
nent sustained by the a 2 -effect in addition to the D-effect. 

The detailed analysis with mode expansion by Legendre 
polynomials on kinetic helicity and magnetic field will be pre¬ 
sented in a subsequent paper in which we focus on the mag¬ 
netic dynamo activity and its dependences on the rotation rate 
and stratification level. Although a quantitative description of 
the dynamo mechanism is beyond the scope of this paper, we 
can speculate from our simulation results that accurate numer¬ 
ical modeling of the kinetic helicity profile in the actual solar 
interior will offer a way to unveil the mystery of the solar 
magnetism. 


5. SUMMARY 

In this paper, the effect of the magnetic field on the mean 
DR profile established in the convection zone was systemati¬ 
cally studied by comparing the MHD and HD models of the 
rotating spherical shell convection in the broad range of the 
initial rotation rate. A fully compressible Yin-Yang MHD 
code which was developed in MYK13 and a stellar model 
consisting of the convection zone and the radiative zone were 
used for the simulation. The critical parameter which controls 
the transition between the anti-solar type and solar-type DR 
profiles was explored with focusing on the “Rossby number 
(Ro) M and the “convective Rossby number (Ro con v)"- In addi¬ 
tion, the bistability of the rotation profile, which has only re¬ 
cently been discovered by Gastine et al. d2014| ), was revisited 
under the influence of the magnetic field. Our main findings 
are summarized as follows. 

1. The critical value of the Ro con v for the transition was 
higher in the MHD model than the HD model. The transi¬ 
tion was more gradual in the model under the influence of the 
magnetic field. Comparing two models at the same Ro con v> 
the magnitude of the DR parameter (a e ) was smaller in the 
MHD model than in the HD model. The magnetic field had a 
crucial role in suppressing the DR and facilitated to produce 
the solar-type rotation profile. Since the Lorentz force of the 
dynamo-generated magnetic field reduced the convective ve¬ 
locity, the MHD model became more rotationally-constrained 
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Fig. 11.— Time and azimuthally-averaged kinetic helicity 0) for the models MR00625, 025, 030, 050, 060, and 150, respectively. The red and blue tones 

denote the positive and negative helicity. The black dotted-line corresponds to the equator. The dashed curve denotes the interface between the convection and 
radiative zones. The upper three models have the anti-solar type DR and the lower three models have the solar-type DR. 
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than the HD model. This would be the reason why the mag¬ 
netic field facilitated to produce the solar-type DR. 

2. The transition from the solar to anti-solar type rotation 
profiles occurred at Ro ~ 1.0 both in the MHD and HD mod¬ 
els. Not only the critical value, the sharpness of the transition 
was also similar between the two models in the view of Ro- 
dependence. The solar-type DR, accompanied by the colum¬ 
nar convection and multi-cell meridional circulation pattern, 
was established in the regime Ro < 1. The anti-solar type 
DR, characterized by the cellular convections and the single¬ 
cell meridional circulation, developed in the regime of Ro > 1 
regardless of the presence of the magnetic field. The transi¬ 
tion between the two rotation profiles was controlled by the 
Ro rather than by the Ro conv . 

3. From the point of view of the Ro conv -dependence, the ro¬ 
tation profile showed, as was observed in earlier studies, the 
bistability near the transition in the HD model while it dis¬ 
apeared when allowing the growth of the dynamo-generated 
magnetic fields except for the model with taking an anti-solar 
type solution as the initial condition. It was found that the 
transition between the anti-solar and solar-type DR profiles 
could be unified in the view of Ro-dependence independently 
of the hysteresis. The transition between the two rotation 
regimes was again better captured by the Ro rather than by 
the Ro conv . 

4. The influence of the magnetic field or the history on the 
convective velocity would be responsible for the dependence 
of the rotation profile on it. Such effects were not reflected 
in the Ro con v by definition. In contrast, it could be definitely 
captured by the Ro. The transition between the solar and anti¬ 
solar type rotation profiles was thus unified in the view of the 
Ro-dependence. 

5. A rich variety of the magnetic dynamo activity could 
be observed in the MHD models. In the regime of the anti¬ 
solar type DR (Ro > 1), the stationary large-scale magnetic 


field with a dipole symmetry was built up both in the stable 
and convection layers except for the model with the smallest 
rotation rate. In contrast, the oscillatory large-scale magnetic 
component was organized in the model with the solar-type 
DR profile (Ro < 1). The latitudinal profile of the kinetic he¬ 
licity would be a key to distinguish the dynamo mechanisms 
between the two rotation regimes. The oscillatory large-scale 
magnetic field could be organized when the strong latitudinal 
gradient of the kinetic helicity existed around the equator. 

Our parametric study was conducted by the simulation 
model with the same value of the Rayleigh number. The su¬ 
percriticality level S (= [Ra-Ra c ]/Ra c ), which is an essential 
measure of the turbulence level, was thus not constant and 
varied with the rotation rate as shown in Table 1 (see the third 
column). In order to compare the models with the supercrit¬ 
icality level fixed, we should change the Ekman number or 
Prandtl number, or density stratification, providing a possible 
difference in the resultant flow pattern. 

Gastine et al. (2013,2014) indicated that, while the transi¬ 
tion between the solar and anti-solar type rotation DR profiles 
would take place at almost the same Ro con v independently of 
the density stratification and Ekman number, the lower Ekman 
number or the stronger density stratification would produce 
the steeper transition. It was also suggested that the higher 
Prandtl number would shifts the critical Ro con v to the higher 
value at least in the high Ekman number regime. Not only on 
the flow properties, these parameters would also have an in- 
fluence on the magnetic dynam o activities (c.f., Kapyla et al. 


|2011| 2014[|Karak et al.|2014[ ). To verify the universality of 
the critical Rossby number for the transition obtained in this 
study, further parameter study is indispensable and is planned 
as our future work. 

The recent development of the astroseismology opens up 
the way to study the large-scale internal flows in the solar- 
type main sequenc e stars with differen t age and thus different 
rotation rate (e.g., [Chaplin et al.|2010 ). Computer simulation 
in tandem with the advanced observation will help deepening 
the understanding of the stellar interior dynamics and stellar 
dynamo activities in the astroseismology era. 
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